Phonon Squeezing in a Superconducting Molecular Transistor 
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Josephson transport through a single molecule or carbon nanotube is considered in the presence of a local 
vibrational mode coupled to the electronic charge. The ground-state solution is obtained exactly in the limit of 
a large superconducting gap, and is extended to the general case by variational analysis. The Josephson current 
induces squeezing of the phonon mode, which is controlled by the superconducting phase difference and by the 
junction asymmetry. Optical probes of non-classical phonon states are briefly discussed. 

PACS numbers: 72.15.Qm. 85.35.Gv. 85.75.-d 



I 

o 



> 
OO 

o 
\o 
o 

-I— > 



I 

o 
o 



X 



Nanoelectromechanical systems (NEMS) offer a way to 
reach the quantum regime of mechanical oscillators. A chal- 
lenging goal consists in creating non-classical vibrational 
states, similar to non-classical states of lightQl with re- 
duced quantum fluctuations in one of the mode quadratures. 
The generation of such states has been suggested in bulk ma- 
terials lisl^. Recently, squeezing scenarios have been sug- 
gested for resonators driven by nonlinear couplings 01 and 
for a Cooper-pair box coupled to a cantilever |6|]. Here we 
consider a superconducting molecular transistor, where the 
molecule carries both electronic and vibrational degrees of 
freedom. We show that squeezing occurs when a Josephson 
cuiTent flows through the molecule, by exploiting the coher- 
ent regime of polaron dynamics [7]. Compared to the fre- 
quency range of micro-cantilevers (less than 1 GHz), here the 
frequency range goes from 1 GHz to 10 THz depending on 
the system (single wall carbon nanotube (SWCNT) or single 
molecule). 

The excitation of molecular vibrations (phonons) by an 
electronic cuiTent has been observed with normal metallic 
leads in several molecules including fullerenes and carbon 
nanotubes The latter are in fact NEMS, where charge 
fluctuations couple to bending |9], stretching |10| or radial 
breathing (REM) modes |11]. Incoherent polaron- like charge 
fluctuations due to transpor t trig ger a non-equilibrium distri- 
bution of the phonon mode lll2ll . In the present paper, we ad- 
dress the superconducting regime. Recently, the proximity ef- 
fect and the Josephson effect 1 13] have been observed in gated 
carbon nanotubes. Yet, in the tunneling regime and for strong 
enough electron-phonon (e-ph) interaction, the Josephson 
current should be suppressed due to the Franck-Condon dress- 
ing factor IT4I1 . In contrast, we show here that for transparent 
lead-molecule contacts, the Josephson effect survives a strong 
e-ph interaction and triggers coherent phonon fluctuations. 
As a striking consequence, the conjugate momentum of the 
molecular distortion displays reduced zero-point fluctuations. 
The superconducting phase difference and the junction asym- 
metry allow to control non-classical squeezed phonon states, 
including nearly (gaussian) minimum-uncertainty states. 

The model Hamiltonian reads H = Hm + Hl + Hr + Ht, 
where Hm, Hlh and Ht respectively describe the uncou- 
pled molecule, two superconducting leads and the molecule- 



lead coupling. Explicitly (we put ?i = fc^ = 1), 



- -) + Un^ni + flb^b, 



Hm^ [eo-A(fe + 5t)] J2 

(1) 

where rio- — d^d^ with da, dj. being the fermion operators 
for spin a =t, J, on the molecular level eo — eo(V^), the po- 
sition of which can be tuned by the gate, U is the repulsive 
Coulomb interaction, and &, 6^ are the boson operators for the 
local mode of frequency 57. The leads are described by stan- 
dard BCS Hamiltonians (j — L,Ris the lead index). 



Hj = J2 i^kCTz + Act,) , ^3k = 



t 



3k A 



(2) 

with the energy dispersion S^k', the Pauli matrices cjx^z act in 
Nambu space. The molecule-lead coupling is given by 



(3) 



where Tj^j^j^ = tja^ g±jo-^0/4^ jg jj^g j^jj lead-molecule 
tunneling amplitude, and (p is the superconducting phase 
difference across the molecule. Averaging over the leads 
yields a partition function for the molecular site, Z[(/)] = 
Tr {e-'3HMyy(^ 0)}, with 



W{(3, 0) = Tr exp ■ 



drdr' (r) S(t - t') (i(r') | 



(4) 

where Tr is the imaginary-time ordering operator, f3 is the 
inverse temperature, E(t — t') = J^j Tj ~ ''"') 

gir) = — J2k i^r + ^kO'z + Acti,)^^ 6{t) is the Green func- 
tion of the uncoupled leads. Assuming a constant normal-state 
density of states in the leads, i^, one obtains: 



s(r-T') = r 



dr 



A I cos - 



Q(t-t') 
(5) 

with Q{t) = 13-^ Eu.^ g-^^„r/^^2 + A2, where uj„ = 
n{2n + l)//3 is a fermionic Matsubara frequency (integer n), 
r = TL + Tr, 7 - (Fl - Tr) /r, and = Tn^e^. 
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First we focus on the case where the relevant molecule dy- 
namics is restricted to the low-frequency domain on the scale 
of the superconducting gap (in particular, F ^ A). In this 
limit, the S term in Eq. becomes local in time with 
Q{t — r') — + A^^(5(t — r') at low frequencies, cj„ ^ A. 
This yields an effective Hamiltonian for the molecule as 



Hef f = Hm + F a' cos — CTj; + 7 sin — CTj 



(6) 



The Hilbert space of i?e// splits into two subspaces, Tii and 
H2, spanned by electron states { | T)) I i) } and {|0d), | Ti)}. 
respectively. In the large A limit, these subspaces are decou- 
pled. We hereafter assume the repulsion U to be relatively 
small, U /2 < F, consistently with good contacts. Therefore, 
the ground state of the system lies in the H2 subspace (3- 
Introducing the Pauli-matrix operators = |0rf)(0rf| — | tJ, 
) (ti I andT+ = (r_)''' = |0d) (ti |, and performing the rota- 
tion Heff = e^^->^/2^e-*^-x/2 with X = arctan (7 tan |^ 
lead to an effective spin-boson Hamiltonian: 



H = nb''b~ [e- X (b + b^)] T^+f{(f>)' 



(7) 



where f = p(0)F, = -y'cos2((/)/2) + 72 sin2(0/2), 
and e — eo + U/2. Any eigenstate of H, and particularly the 
ground state |^o)^ can be written as 



CXj 

^'^/'l^o) = E (^^°^ 10.) +42)1 ti>) In) ,(8) 



ri=0 



where \n) is the n-phonon Fock state. The state (|8} exhibits 
entangled charge and vibrational states. Note that for e = 0, 
the amplitudes in Eq. © fulfill aI°^ = {-ly+^A^nK owing 
to the parity symmetry H , Tx{—1)''^'' = 0. 

The Josephson current flowing through the molecule 
can be obtained as a functional derivative of the gener- 
ating functional = Z[(f) + ^(r)] with respect to the 



source variable ^(r), J 



-2dlnZ^/5mi 



Af- 
ter averaging over the leads and neglecting retardation 
of the kernel E in Eq. (0), the Josephson current re- 
duces to J = Z-iTr|e-^^^^VF(/3,0)j|, where J = 

Frft (— sm((j)/2) (Ta; + 7 cos((?!)/2) ay)dis an effective current 
operator 1 16| in terms of d. In the rotated basis of Eq. Q, the 
current operator takes the form 



J= (F/2p(0)) [(l-72)sin,/)T, + 27Tj,] 



(9) 



Note that in the presence of phonons, the current operator does 
not commute with the Hamiltonian Q. This results in quan- 
tum fluctuations of the current even in the case of e = 7 = 0. 

In the ground state, the relevant experimental quantity is the 
current expectation value, J{(f)) = (^ol^^l^o)' which probes 
the overlap of phonon states in Eq. (|8}. To make further 
progress, we perform a unitary transformation [/„ = e 
on H of Eq. Q, where P = i{b'^ — b) and a = A/fi: 



H' = Vlb'^b - er^ + f [cos(aP) r^^ + sin(aP) Ty] , (10) 



To lowest order in F/J7, for e = 0, one obtains the ground state 
1*0) = e*^/2|a)| Od) - e-'x/2| _ ^^^^^ | ± a) = 

U±a\Qph) are coherent states of the phonon mode. In the limit 
of large a, the Josephson current J{(j)) = Fe~"^| sin((/)/2)| 
becomes strongly suppressed Iil4il . As we show below, in 
a molecular transistor with a relatively strong coupling to 
the leads a sizeable Josephson current may flow and gen- 
erate squeezing without the necessity of an external drive. 
This is in contrast with Ref. j;^ where an external oscilla- 
tor (f2 < IGHz) is coupled to a Cooper-pair box, for which 
fi < F and a < 1 holds. 

In order to probe the generic squeezing properties of the 
molecular phonon, we now compute the mean square fluctu- 
ations of both the displacement coordinate X — b + b'^ and 
the conjugate momentum P in the ground state of the system 
assuming zero temperature. The ground state is found nu- 
merically by truncating the phonon Fock space to a maximum 
number of phonons, Umax ^ A/il. In the range of parameters 
considered, n„iax = 30 is sufficient to ensure the convergence 
of the calculations. When discussing the numerical results, we 
set = 1. 

Coherent charge fluctuations enhance the fluctuations of X, 

6X = ((*o|^^|*o> - (^'o|^|*o>^)^^^ beyond the quan- 
tum zero-point magnitude. Concomitantly, the charge fluc- 
tuations reduce the fluctuations of the momentum, 5P = 
{^o\P^\^oy^^. Note that (1'o|P|^'o) = always, while 
(^'oi^l^'o) = only for e = 0. 

The effective polaronic interaction constant A2/r((/)) of the 
spin-boson Hamiltonian Q depends on the superconducting 
phase difference and becomes infinite at = tt for a symmet- 
ric junction, 7 = 0. Fig. la shows the 0-dependent variation 
of the momentum fluctuation 6P as well as the uncertainty 
SX6P. (The ground-state value of 5X5P for an harmonic os- 
cillator is 1.) Squeezing {5P < 1) occurs for a wide range of 
parameters, and its intensity depends on cf), in accordance with 
the effective polaronic interaction. 

For small F (strong e-ph interaction), the resulting coher- 
ent states I ± a) are not squeezed, i.e., squeezing marks a 
deviation from the displaced oscillator (coherent) states. The 
generation of gaussian squeezing can be understood from the 
Hamiltonian dlOl l by noticing that P^ is not modified by the 
transformation Ua- Expanding the third term in Eq. ilO\ for 
small A gives ^a^T (6^ + {b^)"^) r^, which precisely gener- 
ates two-phonon squeezed states The variation of 6P 
between (p — and tt, where T{(j)) is maximum/minimum, 
shows that squeezing is a crossover phenomenon and is maxi- 
mal for intermediate coupling parameters of the polaron prob- 
lem |3]. In general, squeezing does not involve minimum- 
uncertainty states. Yet, 5X5P ~ 1 can be made very close to 
unity for intermediate A and F (see Fig. lb). 

Squeezed states can only be produced with a sizeable 
Josephson current. Actually, optimal squeezing (in the sense 
of minimum-uncertainty) is generated in the parameter range 
where J is moderately affected by e-ph interaction. To illus- 
trate this, the Josephson current is plotted as a function of the 
bare interaction A^/Fil (Fig. 2a): as expected, J decreases 
but moderately if A is not too large. The inflexion region 
corresponds to the polaron crossover where optimum squeez- 
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FIG. 1: (a) Squeezing of the momentum and (b) Josephson current 
as functions of the phase difference </!>. We take F — 0.5 (circles), 2 
(squares), and 4 (triangles) for A = 0.9, e = 7 = 0. The inset in (b) 
shows the Heisenberg uncertainty as a function of (p. 
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FIG. 2: (a) Josephson current as a function of the electron-phonon 
interaction, at <^ = 7r/2: F — 0.5 (circles), 1 (squares) and 2 (trian- 
gles) for e = 7 = 0; (b) Same as Fig. la, in the asymmetric case: 
7 — 0.5 (squares) and 0.75 (triangles) for A = 1.5, F = 4, e = 0. 



ing is achieved. Interestingly, the asymmetry 7 7^ leads to 
strong (up to 40 percent) and nearly harmonic squeezing, hav- 
ing a weak dependence with the phase difference (Fig. 2b). 

To check that the squeezing property is generic and not only 
restricted to the large A limit, we have treated the case of fi- 
nite A by using a variational ansatz approach |3]. A mod- 
ified polaronic transformation Up{r]) = ^vaib^ -i>)t^ pgj._ 
formed on the original Hamiltonian, Eqs. (l)-(3), followed by 
a harmonic squeezing transformation UgqiP) ~ e^^'''' ^'^ \ 
where 77 and f3 are variational parameters. Averaging over 
the zero-phonon state leads to renormalized tunneling ampli- 
tudes TJ — Tj exp [—^a^rfe^^P^ and the effective interac- 
tion U' ^ U — 277(2 — ?7)(A^/il). Squeezing is measured by 
6P = e~'^^. For large F, using the Hartree approximation, 
one minimizes the ground state energy with respect to r/ and 
f3. This yields a good qualitative agreement with the exact 
large A calculation. For finite A, squeezing is still present 
provided that A > 57. For instance, with the parameters of 
Fig. 2b (A = 1.5, F = 4,7 = 0.75) and = 0, one finds 
SP = 0.79 and 0.76 for A = 100 and 4, respectively. In the 
opposite limit of A < F, assuming A = 0.1, A = 1.2, F = 2, 
7 = 0, one finds (at = tt/2) 6P = 0.77 with SXSP = 1.1, 
indicating that even with rather high frequency modes such as 
RBM in nanotubes ($7 « lOme^), squeezing could be ob- 
tained (assuming Nb contacts, A w ImeV). 

Nonclassical phonon states (other than squeezed states) can 



be generated by projecting the ground state on the cuiTent 
eigenstates, which are in principle accessible experimentally. 
The ground state in the cuiTent state basis (defined from Eq. 

© as J|±) = ±J|±)) reads: |*o) = |*+)|+) + |*-)|-), 
with |$±) = eT»x72^- ^ (e»x/2^W ±g-.x/2^(^2)) |„) 

and x' = arctan ^7 cot In the small F regime, one ob- 
tains, for e = 0, the "phase cat" states, 

|$±) = e^'^'^^ (e*x/2|^^ ^ g-»x/2| _ , (11) 

which are superpositions of opposite coherent states of the 
phonon mode. More generally, states like |$±) are similar to 
those generated in Quantum Optics 1 1]. They caiTy nontrivial 
phases, obtained here by tuning the junction asymmetry 7. In 
particular, for 7 = 0, they are built from phonon Fock states 
with even/odd occupation numbers only, being linked, respec- 
tively, to the |±) current states. At moderate A it is possible to 
achieve such nontrivial states. They have a substantial overlap 
between phonon states and a sizable Josephson current. 

We have also checked that squeezing is robust against en- 
vironmental effects induced by coupling the local vibrational 
mode to an external phonon bath, by using the reduced den- 
sity matrix formalism. Provided that the quality factor Q — 
51/?7, with f] the dissipative coupling constant, is large enough 
(Q ^ 10^ - lO'' has been suggested k^il^ill||]), the squeez- 
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ing property of the ground state is protected due to the gap 
(~ min(ri, A)) in the excitation spectrum of the system. 

A direct probe of squeezed (non-classical) phonon states is 
needed. For low-il bending modes in SWCNT, capacitive de- 
tection through a single-electron transistor can be envisioned. 
Optical detection techniques may be used, such as Resonant 
Raman Scattering (RRS), which has been achieved in carbon 
nanotubes 1 17] for RBM modes. In RRS, the incident pho- 
ton excites an electronic transition within the molecule, and a 
photon is re-emitted with excitation of the phonon modes. We 
denote d! the state corresponding to a low-lying molecular or- 
bital, which is assumed to be decoupled from the molecular 
vibrations. The molecule Hamiltonian then reads: 

H'm = Hm + t'd'W + {Ce'^^^dW + H.c) . (12) 

Using the golden rule to calculate the Raman transition 
rate, we assume that the initial state is the projection 
of the ground state (|8j on the zero-electron subspace, 
\i) = |$(o))|0)d| Ti)rf„ with the phonon part |$(°)> = 
^^P^i^'ln), while the final state with one electron in 
the orbital d (together with a hole in d'), involves an arbi- 
trary number of phonons without a net displacement, |/) = 
|'^)|f)ci| ~ cr)d'. Resonant Raman emission lines (Stokes) 
will appear at energies eo — e' + nil — Eo{(f>) (EQ{(j>) is the 



ground state energy), with probabilities given by (n|$(°^)^ = 

I Ar?"* p; this yields a spectroscopy of the squeezed state |<I>(°)). 
While absorption from a filled molecular orbital {d') probes 
1$^*'^), similarly absorption towards an empty orbital will 

probe l^'^'^'f) = J2'^=o ^"^^ I") taking out an electron from 
the state | ti)d- Concerning the "cat states", experimental 
evidence could be gained in principle by far-infrared optical 
absorption involving a transition from the subgap (Andreev) 
bound states to quasiparticle states in the leads. 

In summary, we have studied theoretically the generation of 
phonon squeezing by coherent charge fluctuations in a super- 
conducting molecular junction. Squeezing occurs for a wide 
range of parameters and is maximal in the polaron crossover 
regime. A nearly minimum-uncertainty state with about 40 
per cent squeezing can be obtained provided that the local vi- 
brational mode is weakly coupled to the environment. For 
these purposes, suspended nanotubes present the advantage 
that their various vibrational modes are well characterized and 
the corresponding quality factors are large. Finally, we have 
discussed the possibility of using optical techniques to detect 
the spectral properties of non-classical phonon states. 
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